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1.  INTRODUCTION 


In  the  study  of  nonlinear  acoustic  propagation  through  an  inhomogeneous  fluid 
such  as  seawater,  it  is  helpful  to  have  explicit  relations  for  the  acoustic  properties 
of  the  fluid  as  a  function  of  temperature,  pressure,  and  composition.  For  both  fresh 
water  and  seawater,  several  empirical  relations  give  the  sound  speed  c  and  the  density 
p  as  functions  of  temperature,  pressure,  and,  in  the  case  of  seawater,  salinity  S. 
We  are,  however,  unaware  of  any  similar  relations  for  the  coefficient  of  nonlinearity 
/ 3 .  Endo  [1]  has  tabulated  /?  for  seawater  over  a  broad  range  of  temperatures  and 
pressures  for  salinities  of  25,  30,  35,  and  40  parts  per  thousand  (ppt).  For  computer 
implementations  of  nonlinear  propagation  through  seawater,  however,  relations  with 
explicit  temperature,  pressure,  and  salinity  dependence  are  preferable  over  tables  of 
values.  Endo’s  tabular  values  were  useful  in  checking  our  calculations. 

In  this  work  three  relations  for  /?  are  developed:  one  for  fresh  water  as  a  function 
of  temperature  and  pressure,  one  for  seawater  with  a  salinity  of  35  ppt  as  a  function  of 
temperature  and  pressure,  and  one  for  seawater  as  a  function  of  temperature,  pressure, 
and  salinity.  Also  developed  is  a  similar  set  of  relations  for  the  term  A  =  0/(pc5)1^2. 

The  reason  for  developing  separate  relations  for  A  is  that  it  arises  as  an  important 
variable  in  nonlinear  geometrical  acoustics  [2,3].  Under  the  assumptions  of  nonlinear 
geometrical  acoustics,  the  equation  governing  the  propagation  of  a  progressive  finite- 
amplitude  acoustic  wave  in  an  inhomogeneous  fluid  may  be  converted  to  the  equation 
governing  the  propagation  of  a  similar  plane  wave  in  a  homogeneous  fluid.  The 
conversions  required  are  the  following: 


where  p  is  the  acoustic  pressure,  x  and  p  are.  respectively,  the  equivalent  plane  wave 
propagation  distance  and  pressure.  A  is  the  ray  tube  area,  and  (  is  the  ray  path  length. 
The  subscript  0  indicates  conditions  at  the  source  point.  Equation  (1.2)  makes  clear 
why  an  explicit  relation  for  A  as  a  function  of  temperature,  pressure,  and  salinity  is 
helpful. 
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2.  OVERVIEW  OF  PROCEDURE 


The  coefficient  of  nonlinearity  is  defined  as  follows:  /?  =  B/2A  + 1,  where  A  and  B 
are  the  first  and  second  coefficients  of  the  Taylor  series  expansion  of  the  total  pressure 
P  in  terms  of  the  relative  density  change.  (In  fresh  water  at  room  temperature  and 
atmospheric  pressure,  /?  has  a  value  of  approximately  3.5.)  Rudnick  [4]  has  shown 
that  B/ A  may  be  expressed  in  terms  of  quantities  that  are  easily  related  to  physical 
data.  Use  of  the  expression  for  B/ A  as  shown  in  the  paper  by  Beyer  [5]  gives 


P 


(2.1) 


where  the  subscripts  indicate  constant  conditions,  and  T,  a,  and  Cp  are,  respectively, 
the  absolute  temperature,  the  coefficient  of  thermal  expansion,  and  the  specific  heat  at 
constant  pressure.  Note  that  the  first  and  second  terms  in  Eq.  (2.1)  contain  different 
partial  derivatives  of  the  sound  speed  c(P,  T)  evaluated  at  constant  composition.  The 
term  that  contains  (dc/dP)p  is  hereinafter  referred  to  as  “term  1”,  whereas  the  term 
that  contains  (dc/OT)p  is  hereinafter  referred  to  as  “term  2".  Rudnick  noted  that, 
while  all  parts  make  a  significant  contribution  to  /3,  term  1  is  numerically  larger  than 
term  2  in  all  common  fluids. 


In  our  f  valuation  of  [i .  published  empirical  relations  were  used  as  much  as  possible. 
When  a  derivative  of  an  acoustical  property  was  required,  the  derivative  was  obtained 
by  differentiating  the  published  empirical  relation  analytically,  then  evaluating  the 
result  numerically. 

To  the  best  of  our  knowledge,  there  is  only  one  set  of  empirical  relations  for  the 
acoustical  properties  of  fresh  water  that  is  valid  over  a  wide  range  of  temperatures 
(0  1 00°  ( ' )  and  pressures  (1  lOllO  bars  absolute).  The  relations  used  for  the  sound 
speed  and  density  ol  fresh  water  arc’  from,  respectively,  Chen  and  Millero  [(»]  and 
Chen  <  t  U.I.  [7],  The  estimate’s  of  the*  coefficient,  of  thermal  expansion  were  obtained 
using  p  from  Chen  cl  til.  [7]  in  t  lie  expression  —  f  lor  Cp.  on  the  other  hand, 

we  used  a  polynomial  lit  to  a  liesli  water  (  '/•  data  set  (8). 

I  he  situation  is  slightly  dillerent  lor  seawater:  at  least  two  published  empirical 
relations  that  are  valid  over  a  wide  range  of  temperatures  (0  1()°C)  and  pressures 
(1  800  bars  absolute)  exist  for  both  the  sound  speed  and  the’  density.  For  the1  sound 
speed  the  choice  is  between  the  relation  developed  by  Lovett  [!)]  and  that  developed 
by  Chen  and  Millero  1 1  or  For  the’  density  the’  choice’  is  between  the  expressions 


developed  by  Gebhart  and  Mollendorf  [11]  and  those  from  Chen  and  Millero  [12]. 
The  choice  for  the  coefficient  of  thermal  expansion  is  somewhat  dependent  on  the 
choice  made  for  the  density.  The  choice  is  either  the  evaluation  of  —  ^  (ff)p5  using 

p  from  Gebhart  and  Mollendorf  [11]  or  the  explicit  relation  for  a  given  by  Chen  and 
Millero  [12].  Our  estimate  of  the  specific  heat  at  constant  pressure  in  seawater  Cp  was 
obtained  by  summing  the  estimate  of  Cp  at  one  atmosphere,  which  was  calculated 
using  a  relation  from  Millero  et  al.  [13],  with  the  contribution  from  the  isothermal 
integration  of  [14] 


The  estimates  of  a  and  p  for  Eq.  (2.2)  were  obtained  from  the  appropriate  aforemen¬ 
tioned  references. 


3.  UNCERTAINTY  IN  ESTIMATES  OF  /?  IN  SEAWATER 


To  help  decide  which  empirical  relations  to  use  for  c,  p,  and  a  in  seawater  and, 
more  importantly,  to  estimate  the  uncertainty  in  the  calculation  of  /?,  we  numerically 
implemented  all  of  the  aforementioned  empirical  relations.  The  absolute  difference  in 
(3  due  to  using  two  different  empirical  relations  was  calculated  for  one  property  at  a 
time.  The  difference  in  /3  caused  by  terms  1  and  2  was  calculated  separately  to  see 
which  term  in  /3  is  more  sensitive. 

Term  1,  the  term  in  f3  that  contains  (dc/dP) j,  was  found  to  be  the  largest  source 
of  uncertainty.  Estimates  of  the  uncertainty  were  obtained  by  calculating  term  1  using 
Lovett’s  relation  for  c  [9]  while  holding  all  other  relations  constant.  (For  p  and  a,  the 
previously  referenced  Chen  and  Millero  relations  were  used.)  This  result  was  then 
subtracted  from  a  calculation  of  term  1  obtained  using  Chen  and  Millero’s  relation 
for  c  [10].  The  difference  was  calculated  over  a  realistic  range  of  temperatures  and 
pressures,  that  is,  a  range  of  temperatures  and  pressures  that  might  be  found  in  an 
ocean:  a  temperature  range  of  —4-40°  for  pressures  of  1  and  100  bars  absolute,  but 
a  temperature  range  of  only  —4-10°  for  pressures  of  200,  300,  . . .,  800  bars  absolute. 
Shown  in  Fig.  3.1  is  a  plot  of  the  difference  in  term  1  for  salinities  of  30  and  40  ppt. 
(Solid  lines  are  for  salinity  of  30  ppt,  dashed  for  40  ppt.)  For  these  salinities,  the 
difference  is  always  less  than  0.1,  typically  0.05.  For  a  salinity  of  35  ppt  (not  shown), 
the  absolute  difference  is  always  less  than  0.05  and  typically  less  than  0.025.  For 
contrast,  Fig.  3.2  shows  a  plot  of  the  corresponding  difference  in  term  2,  which  was 
calculated  in  an  identical  fashion.  The  magnitude  of  these  differences  is  much  smaller, 
typically  less  than  0.0025. 

The  uncertainties  in  f3  that  stem  from  the  different  empirical  relations  for  p  and  a 
were  examined  in  the  same  manner  as  those  for  c..  Both  terms  1  and  2  were  evaluated 
twice,  once  using  the  density  relation  from  Gebhart  and  Mollendorf  [11]  and  once 
using  Chen  and  Millero’s  relation  [12],  Because  o  is  closely  related  to  p ,  it  was  varied 
simultaneously.  (In  all  cases,  the  relation  used  for  c  was  that  of  Chen  and  Millero 
[10].)  Plots  of  the  differences  generated  by  terms  1  and  2  are  shown,  respectively, 
in  Figs.  3.3  and  3.4.  The  largest  differences  in  ,i  that  result  are  in  term  2,  and  the 
differences  are  always  less  than  0.01. 

Finally,  a  check  on  the  pressure  correction  term  in  C'p,  Fq.  (2.2),  was  made  at 
zero  salinity.  Between  1  and  1000  bar,  the  pressure  increment  in  C/>  given  by  Fq.  (2.2) 
deviated  by  around  10%  from  the  values  in  Ref  S.  Nevertheless,  the  corresponding 
nncertainity  in  >3.  through  term  2,  is  only  1  part  in  lOF 
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FIGURE  3.1 

PLOTS  OF  THE  ABSOLUTE  DIFFERENCE  IN 
(J  versus  TEMPERATURE  FOR  ABSOLUTE  PRESSURES 
1, 100,  200, 300  .800  BAR  CAUSED  BY  USING  DIFFERENT 
EMPIRICAL  RELATIONS  FOR  c  IN  TERM  1 


FIGURE  3.2 

PLOTS  OF  THE  ABSOLUTE  DIFFERENCE  IN 
3  versus  TEMPERATURE  FOR  ABSOLUTE  PRESSURES 
OF  1, 100, 200,  300  . .800  BAR  CAUSED  BY  USING  DIFFERENT  • 

EMPIRICAL  RELATIONS  FOR  C  IN  TERM  2 
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FIGURE  3.3 

PLOTS  OF  THE  ABSOLUTE  DIFFERENCE  IN 
3  versus  TEMPERATURE  FOR  ABSOLUTE  PRESSURES 
OF  1, 100,  200, 300.  .800  BAR  CAUSED  BY  USING  DIFFERENT 
EMPIRICAL  RELATIONS  FOR  p  AND  a  IN  TERM  1 


FIGURE  3  4 

PLOTS  OF  THE  ABSOLUTE  DIFFERENCE  IN 
3  versus  TEMPERATURE  FOR  ABSOLUTE  PRESSURES 
1,  100,  200,  300.  800  BAR  CAUSED  BY  USING  DIFFERENT 
EMPIRICAL  RELATIONS  FOR  p  AND  a  IN  TERM  2  • 
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Rather  than  embark  on  an  investigation  of  the  origin  of  the  above  cited  differences, 
we  simply  chose  to  use  the  relations  from  Chen,  Millero,  and  their  co-workers  over 
those  of  other  workers.  Our  reasoning  was  uniformity  of  approach:  Chen,  Millero, 
and  their  co-workers  have  developed  empirical  relations  not  only  for  c,  but  also  for 
p,  a,  and  Cp.  Thus,  for  seawater,  the  relations  in  Refs.  10  and  12  were  used  for  c 
and  p.  Estimates  of  a  were  obtained  using  the  explicit  relation  for  a  in  seawater 
from  Ref.  12.  When  using  our  polynomial  expressions  to  calculate  f3,  however,  one 
must  always  bear  in  mind  that  the  differences  cited  above  are,  generally  speaking, 
the  accuracy  limits  in  our  results. 
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4.  RESULTS 


4.1  PLOTS  OF  c,  p,  ft,  AND  A 

Plots  of  c,  p,  ft,  and  A  for  fresh  water  are  shown  in  Figs.  4. 1-4.4  and  for  seawater 
in  Figs.  4. 5-4. 8.  For  fresh  water,  the  acoustic  properties  are  plotted  as  functions  of 
temperature  for  the  absolute  pressures  of  1,  100,  200,  300,  ...,  1000  bars.  Because 
the  fresh  water  Cp  relation  is  valid  up  to  only  90°C,  our  relations  for  (3  and  A  in  fresh 
water  must  be  used  with  caution  above  90°C.  The  curves  are  therefore  dashed  in  this 
region.  In  the  case  of  seawater,  Figs.  4. 5-4. 8,  the  salinity  is  held  constant  at  35  ppt. 
The  acoustic  properties  are  plotted  as  functions  of  temperature  for  absolute  pressures 
of  1,  100,  200,  300,  . . .,  800  bars.  Because  the  accuracy  of  the  relations  for  c  and  p 
in  seawater  is  uncertain  below  0°C,  the  values  for  c  and  p  below  0°C  are  dashed  lines 
in  this  region.  Similarly,  our  derived  relations  for  ft  and  A  in  seawater  must  be  used 
with  caution  below  0°C,  and  they  too  are  dashed  lines  in  this  region. 

Some  discussion  of  the  plots  is  in  order.  The  plots  of  the  sound  speed  and  density 
for  both  fresh  water  and  seawater  are  as  expected:  the  density  decreases  with  rising 
temperature,  but  increases  with  rising  pressure,  whereas  the  sound  speed  increases 
with  pressure  and,  up  to  around  70-80°,  increases  with  temperature  also.  The  plots  of 
the  coefficient  of  nonlinearity  reveal  a  more  complex  temperature  and  pressure  depen¬ 
dence.  In  seawater,  (3  increases  with  pressure  for  the  temperature  range  considered. 
However,  in  fresh  water,  the  sign  of  the  pressure  derivative  of  (3  is  temperature  de¬ 
pendent.  Below  about  60°C,  (3  increases  monotonically  with  pressure,  whereas  above 
that  temperature,  ft  initially  decreases  with  rising  pressure.  For  both  fresh  water  and 
seawater,  the  pressure  effect  on  ft  tends  to  decrease  with  rising  temperature.  The 
variation  in  ft  near  0°C  is  about  25%,  whereas  at  higher  temperatures  the  variation 
decreases  to  about  5%.  In  contrast  with  ft,  the  pressure  effect  on  A  increases  with 
temperature  for  both  fresh  water  and  seawater.  Near  0°C,  the  variation  is  about  5%, 
but  by  40°C  the  variation  has  increased  to  about  25%,  and  (in  fresh  water)  by  90°(' 
the  variation  has  increased  to  about  40%.  Aside  from  a  deviation  at  low  pressure 
and  temperature,  the  temperature  and  pressure  dependences  of  A  for  both  seawater 
and  fresh  water  are  generally  opposite  those  of  the  sound  speed  r.  Phis  is  reasonable 
given  that  ch^2  appears  in  t he  denominator  ol  A. 

A  comparison  of  the  isobar  plots  of  ft  for  fresh  water  and  seawater  (Fig.  4.3  with 
Fig.  1.7)  casts  some  doubt  on  the  validity  of  our  extrapolation  of  the  underlying 
expressions  for  the  acoustic  properties  below  0°('.  Note  that,  in  our  fresh  water  ft 
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FIGURE  4.1 

PLOTS  OF  SOUND  SPEED  versus  TEMPERATURE 
FOR  VARIOUS  ABSOLUTE  PRESSURES  (BAR) 
IN  FRESH  WATER 


FIGURE  4.2 

PLOTS  OF  DENSITY  versus  TEMPERATURE 
FOR  VARIOUS  ABSOLUTE  PRESSURES  (BAR) 
IN  FRESH  WATER 


AS-91  -653 


AS-91  <54 


FIGURE  4.3 

PLOTS  OF  3  versus  TEMPERATURE 
FOR  VARIOUS  ABSOLUTE  PRESSURES  (BAR) 
IN  FRESHWATER 


FIGURE  4  4 

PLOTS  OF  A  versus  TEMPERATURE 
FOR  VARIOUS  ABSOLUTE  PRESSURES  (BAR) 
IN  FRESHWATER 


AS  91  -655 


AS  91  656 
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FIGURE  4.5 

PLOTS  OF  SOUND  SPEED  versus  TEMPERATURE 
FOR  VARIOUS  ABSOLUTE  PRESSURES  (BAR) 

IN  SEAWATER  WITH  A  SALINITY  OF  35  ppt  [1 0] 

AS-91-657 


TEMPERATURE  -  °C 


FIGURE  4.7 

PLOTS  OF  p  versus  TEMPERATURE 
FOR  VARIOUS  ABSOLUTE  PRESSURES  (BAR) 
IN  SEAWATER  WITH  A  SALINITY  OF  35  ppt 


FIGURE  4.6 

PLOTS  OF  DENSITY  versus  TEMPERATURE 
FOR  VARIOUS  ABSOLUTE  PRESSURES  (BAR) 

IN  SEAWATER  WITH  A  SALINITY  OF  35  ppt  [1 2] 

AS-91-65* 


FIGURE  4.8 

PLOTS  OF  A  versus  TEMPERATURE 
FOR  VARIOUS  ABSOLUTE  PRESSURES  (BAR) 
IN  SEAWATER  WITH  A  SALINITY  OF  35  ppt 


AS-91  659 
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results,  the  /?  curves  diverge  rapidly  as  they  approach  the  freezing  point,  whereas  in 
our  seawater  (3  results,  the  f3  curves  do  not  exhibit  a  similar  divergence  below  0°C. 
The  reason  may  be  our  unjustified  extrapolation  of  the  underlying  expressions  for  c, 
p,  a,  and  Cp  in  seawater.  Although  these  relations  are  not  stated  as  being  accurate 
below  0°C,  they  are  the  best  available.  Because  we  believe  that  having  a  expression 
for  /?  with  an  uncertain  accuracy  below  0°C  is  preferable  over  no  expression  at  all, 
we  extrapolated  the  underlying  expressions  to  develop  our  relation  for  (3. 


4.2  POLYNOMIAL  EXPRESSIONS  FOR  f3  AND  A 


Our  polynomial  expressions  for  (3  and  A  are  now  presented.  The  order  of  our  poly¬ 
nomial  fit  was  chosen  to  be  similar  to  that  used  in  the  empirical  fit  to  the  dominant 
underlying  acoustic  property,  the  sound  speed  c.  The  sound  speed  relation  given  by 
Chen  and  Millero  [7]  is  fifth  order  in  temperature,  third  order  in  pressure,  and  second 
order  in  salinity.  The  order  of  the  salinity  terms  was,  however,  reduced  from  second 
to  first  order  because  the  sound  speed  relation  given  by  Chen  and  Millero  [7]  is  valid 
for  a  broader  range  of  salinities  than  is  required  (5-40  ppt  rather  than  our  30-40  ppt) . 
The  resulting  fit  of  the  polynomials  to  the  original  data  is  more  than  adequate.  The 
quality  of  the  fits  of  the  polynomials  to  the  original  data  was  assessed  through  direct 
comparison,  and  the  comparison  was  performed  over  the  claimed  ranges  of  tempera¬ 
tures,  pressures,  and  salinities.  The  difference  between  the  original  data  and  the  fit 
is  always  far  less  than  the  previously  discussed  uncertainty  in  the  overall  calculation. 

Our  polynomial  expressions  for  f3  were  obtained  as  follows.  For  seawater,  Eq.  (2.1 ) 
was  evaluated  over  the  following  range:  temperature  —4  to  40°C,  pressure  1-800  bars 
absolute,  and  salinity  30-40  ppt.  The  resulting  data  points  were  then  fitted  to  a 
polynomial  expression  that  is  fifth  order  in  temperature,  third  order  in  pressure,  and 
linear  in  salinity.  The  polynomial  expression  is 

P=  E  Blmn(S-M)‘P'nTn  ,  (4.1) 

/,m,n 

where  the  salinity  S  is  in  ppt,  the  temperature  T  is  in  °C,  and  the  pressure  P  is  bars 
gauge,  that  is,  relative  to  1  atm  =  1.01325  bar.1  The  coefficients  Bimn  are  listed  in 
Table  4.1.  The  maximum  absolute  error  in  the  fit  is  0.0015,  and  the  typical  error  is 
about  half  that.  (Check  value:  For  a  salinity  of  32.5  ppt ,  a  pressure  of  100  bar  gauge, 
and  a  temperature  of  J5°C,  !i  is  3.672556.) 


The  act  ual  salinity  of  t  he  ocean  water  column  is  frequently  unknown.  In  its  place, 
people  often  assume  a  salinity  ol  35  ppt  because  this  salinity  is  so  common  in  the 

'This  is  tile  convention  used  by  most  authors,  including  Chen,  Millero,  and  their  co-workers,  for 
presenting  empirical  fits  to  seawater  and  fresh  water  property  data. 


TABLE  4.1 

Coefficients  in  polynomial  expression  for 


l 

m 

n 

$lmn 

“d3 5 - 

^mn 

0 

0 

0 

+3.122127 

+3.44166 

+3.44223 

0 

0 

1 

+0.0262768 

+0.0176224 

+0.0176264 

0 

0 

2 

-0.000439164 

-0.000478735 

-0.00047893 

0 

0 

3 

+0.00000520032 

+0.0000129976 

+0.0000130006 

0 

0 

4 

-3.07095  x  10-8 

-1.85398  x  10"7 

-1.85429  x  10-7 

0 

0 

5 

+6.12344  x  10~u 

+3.36832  x  10“ 10 

+3.37367  x  lO"10 

0 

1 

0 

+0.00139138 

+0.00104125 

+0.0010417 

0 

1 

1 

-0.000051 12S 

-0.0000457645 

-0.0000457763 

0 

1 

2 

+7.40662  x  10“7 

+5.81691  x  10-7 

+5.85248  x  10-7 

0 

1 

3 

-6.41226  x  10~y 

+5.16276  x  10-9 

+5.14004  x  lO"9 

0 

1 

4 

+2.53274  x  10- 11 

—3.81942  x  10-" 

-3.82208  x  10- 11 

0 

2 

0 

-2.15517  x  10-7 

-1.69183  x  10“7 

-1.69416  x  lO"7 

0 

2 

1 

+3.45228  x  10-9 

+  1.69021  x  10"8 

+  1.69053  x  lO"8 

0 

2 

2 

+  1.44653  x  10~10 

-3.37325  x  10~10 

-3.37744  x  lO"10 

0 

2 

3 

-2.82686  x  10- 12 

-3.82096  x  10- 12 

-3.79902  x  10~12 

0 

2 

4 

+  1.26217  x  10-" 

+4.52996  x  10-“ 

+4.49897  x  10-“ 

0 

3 

0 

-1.12614  x  10-'° 

-6.58146  x  10-" 

-6.57756  x  10" 11 

0 

3 

1 

+4.51141  x  10-12 

+6.60716  x  10-'2 

+6.60637  x  10-12 

0 

1 

1 

1 

1 

1 

1 

1 

i 

1 

3 

0 

0 

0 

1 

1 

1 

> 

2 

2 

2 

0 

1 

2 

0 

1 

2 

0 

1 

•  > 

-2.91434  x  10-“ 

-1.63621  x  10-“ 
+0.00833572 
-0.000232107 
+0.00000107536 
-0.00001 15765 
+  7.2387  x  1()“9 

4  1  0786  v  1()-8 
+  3.2689  x  H)-9 
41.6 156  .  10  10 

1 .807  1 1  ■  10  " 

-1.63501  x  lO" 13 

world’s  oceans.  For  these  reasons,  a  less  complicated  relation  for  (3  for  the  salinity  of 
35  ppt  was  also  developed, 


/J35  =  Y  B3JnPmTn 

m,n 


(4.2) 


Valid  over  the  same  range  of  temperatures  and  pressures,  its  coefficients  are  denoted 
Bmn  and  are  listed  in  Table  4.1.  The  maximum  absolute  error  in  the  fit  is  0.0001. 
(Check  value:  For  a  pressure  of  98.98675  bar  gauge  and  a  temperature  of  15°C,  /?35 
is  3.683292.) 

For  fresh  water,  Eq.  (2.1)  was  evaluated  over  the  temperature  range  0-100°C 
and  the  pressure  range  1-1000  bars  absolute.  A  polynomial  fit  of  the  same  order  in 
temperature  and  pressure  as  above  was  then  obtained, 

/?°  =  Y  B°mnpmTn  •  (4-3) 

m  ,n 


The  coefficients  B^n  are  also  listed  in  Table  4.1,  and  the  maximum  absolute  error  of 
the  fit  is  0.0005.  (Check  value:  For  a  pressure  of  98.98675  bar  gauge  and  a  temperature 
of  15°C,  /9°  is  3.508381.) 

Polynomial  fits  for  A  =  /3/(pc5)1^  that  are  valid  over  the  same  ranges  of  tem¬ 
peratures  and  pressures  for  seawater  with  salinities  of  30-40  ppt,  seawater  with  a 
salinity  of  35  ppt,  and  fresh  water  were  developed  as  above.  The  relations  are  shown 
in  Eqs.  (4.4)-(4.6),  and  the  coefficients,  which  are  denoted,  respectively,  Z./mn,  /'3^„, 
and  L°mn,  are  listed  in  Table  4.3: 


A  =  [Y  L‘mn(S  -35)lPmTn\  X  10-9  s5/2m-'kg-1/2  ,  (4.4) 

\t.m,n  ) 


l.r>  _  1 

[ZL:*nPm'r) 

\m,n  / 

X  10~9 

s^m-'kg-'/2  , 

(4.5) 

0  —  1 

[S  C 

\m,u  / 

x  10“9 

s5/2m-'kg-,/2  . 

(4.6) 

The  maximum  absolute  errors  in  these  fits  are,  respectively,  0.0004  x  10~9,  0.0001  x 
10-9,  and  0.0015  x  10-9.  (Check  values:  For  a  salinity  of  32.5  ppt,  a  pressure  of 
100  bar  gauge,  and  a  temperature  of  15°C,  A  is  1.270685  x  10~9.  For  a  pressure 
of  98.98675  bar  gauge  and  a  temperature  of  15°C,  A35  is  1.267513  x  10-9  and  A°  is 
1.309151  x  10"9.) 
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TABLE  4.2 

Coefficients  in  polynomial  expression  for  A. 


/ 

m 

n 

T3 

mn 

Limn 

- 

^mn 

0 

0 

0 

+  1.340649 

+  1.34255 

+  1.34282 

0 

0 

1 

-0.00123676 

-0.00369707 

-0.00370087 

0 

0 

2 

-0.0000219506 

-0.0000454961 

-0.0000454702 

0 

0 

3 

+7.51912  x  10"7 

+0.00000411388 

+0.0000041155  | 

0 

0 

4 

-2.82344  x  10“9 

-7.30211  x  10“8 

-7.30791  x  10“8 

0 

0 

5 

-5.22167  x  10"12 

+2.49154  x  10-10 

+2.49717  x  10~10 

0 

1 

0 

+0.000197704 

-2.88818  x  10-7 

-4.62069  x  10"7 

0 

1 

1 

-0.0000256345 

-0.0000191962 

-0.0000191921 

0 

1 

2 

+4.70095  x  10~7 

+4.28444  x  10~7 

+4.28422  x  lO"7 

0 

1 

3 

-4.83161  x  10-9 

-4.12041  x  10~9 

-4.1215  x  lO'9 

0 

1 

4 

+  1.90758  x  10-11 

+4.36867  x  10-11 

+4.37169  x  lO-11 

0 

2 

0 

-2.69563  x  10"7 

-1.43851  x  10-7 

-1.43841  x  lO"7 

0 

2 

l 

+  1.46666  x  10-8 

+  1.65626  x  10~8 

+  1.65666  x  10-8 

0 

2 

2 

-1.99597  x  10-10 

-3.94804  x  lO-10 

-3.95058  x  lO" 10 

0 

2 

3 

+  1.59563  x  10"12 

+2.14739  x  10~12 

+2.15271  x  10"12 

0 

2 

4 

-5.86015  x  10"15 

-1.49805  x  10'' 14 

-1.50419  x  10" 14 

0 

3 

0 

+4.01117  x  10-11 

+2.48433  x  10"u 

+2.48674  x  10~ 11 

0 

3 

1 

-1.60177  x  10-12 

-2.5426  x  10~12 

-2.54489  x  10"12 

0 

3 

2 

+8.41646  x  10~15 

+4.69605  x  10~14 

+4.70247  x  10-14 

1 

0 

0 

-0.000374879 

1 

0 

1 

-0.0000614364 

1 

0 

2 

+4.68875  x  10-7 

1 

1 

0 

-0.00000518082 

1 

1 

1 

+  1.53554  x  10~7 

1 

1 

2 

+9.56226  x  lO"10 

1 

2 

0 

+2.92378  x  10~9 

1 

2 

1 

+2.13857  x  10_M 

1 

2 

2 

-3.87438  x  10“ 12 

4.3  DISCREPANCY  WITH  THE  TABULATED  VALUES  OF  ENDO  [1] 


A  small  discrepancy  exists  between  our  results  and  those  of  Endo  [lj.  The  magni¬ 
tude  of  the  difference  in  terms  of  j3  is  at  most  0.15  at  high  temperature  and  pressure. 
The  results  for  high  pressure,  low  temperature  as  well  as  those  for  low  pressure,  high 
temperature  appear  to  agree  well.  Thus  a  cross  term  appears  to  be  the  source  of  the 
problem.  However,  despite  numerous  tests  and  repeated  checking,  the  source  of  the 
discrepancy  could  not  traced.  The  checks  included  evaluating  the  sound  speed 
derivatives,  Eq.  (2.1),  numerically  for  comparison  with  the  values  obtained  via  ana¬ 
lytical  differentiation.  The  agreement  was  within  1  ppt  for  term  1,  and  even  less  for 
term  2. 

The  discrepancy  is  particularly  unsettling,  as  both  our  and  Endo’s  underlying 
expressions  for  c,  p,  and  a  are  the  same.  However,  the  discrepancy  is  small,  of  the 
same  order  as  the  uncertainty  in  0  caused  by  using  the  different  expressions  for  c,  p, 
and  a. 


5.  SUMMARY 


In  this  work  three  polynomial  relations  for  0  are  developed:  one  for  fresh  water 
as  a  function  of  temperature  and  pressure,  one  for  seawater  with  a  salinity  of  35  ppt 
as  a  function  of  temperature  and  pressure,  and  one  for  seawater  as  a  function  of 
temperature,  pressure,  and  salinity.  Also  developed  are  a  similar  set  of  relations  for 
the  term  A  =  0/(pc5y/2.  Plots  of  the  sound  speed  c,  the  density  p,  as  well  as  0  and 
A  for  both  fresh  water  and  seawater  with  a  salinity  of  35  ppt  are  presented  for  a  wide 
range  of  temperatures  and  pressures.  Although  the  precision  of  the  polynomial  fits 
to  the  underlying  data  is  quite  good  (typically  better  than  1  ppt),  the  use  of  different 
underlying  property  data  leads  to  uncertainities  in  0  of  about  ±3%.  Due  to  a  lack  of 
comparison  information,  a  similar  estimate  of  the  uncertainty  is  not  available  for  0 
in  fresh  water. 
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